LAPLACIANS ON THE BASILICA JULIA SET 



LUKE G. ROGERS AND ALEXANDER TEPLYAEV 

Abstract. We consider the basilica Julia set of the polynomial P(z) = 
z 2 — 1 and construct all possible resistance (Dirichlet) forms, and the cor- 
responding Laplacians, for which the topology in the effective resistance 
metric coincides with the usual topology. Then we concentrate on two 
particular cases. One is a self-similar harmonic structure, for which the 
energy renormalization factor is 2, the spectral dimension is log 9/ log 6, 
and we can compute all the eigenvalues and eigenfunctions by a spectral 
decimation method. The other is graph-directed self-similar under the map 
z i— > P(z); it has energy renormalization factor y/2 and spectral dimension 
4/3, but the exact computation of the spectrum is difficult. The latter 
Dirichlet form and Laplacian are in a sense conformally invariant on the 
basilica Julia set. 
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1. Introduction 

In the rapidly developing theory of analysis on fractals, the principal ex- 
amples are finitely ramified self-similar fractal sets that arise as fixed points 
of iterated function systems (IFS). For example, the recent book of Strichartz 
[22] gives a detailed account of the rich structure that has been developed 
for studying differential equations on the the well known Sierpinski Gasket 
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fractal, primarily by using the methods of Kigami (see [H]). Some general- 
izations exist to fractals generated by graph-directed IFS and certain random 
IFS constructions [3 O E], but it is desirable to extend the methods to other 
interesting cases. Among the most important and rich collections of fractals 
are the Julia sets of complex dynamical systems (see, for instance, [31 H| [16]). 
In this paper we construct Dirichlet forms and Laplacians on the Julia set 
of the quadratic polynomial P(z) = z 2 — 1, which is often referred to as the 
basilica Julia set, see Figured! The basilica Julia set is particularly interesting 
because it is one of the simplest examples of a Julia set with nontrivial topol- 
ogy, and analyzing it in detail shows how to transfer the differential equation 
methods of |22j to more general Julia sets. 




Figure 1. The basilica Julia set, the Julia set of z 2 — 1. 

Another reason for our interest in the basilica Julia set comes from its ap- 
pearance as the limit set of the so-called basilica self-similar group. This class 
of groups came to prominence because of their relation to finite automata and 
groups of intermediate growth, first discovered by Grigorchuk. The reader can 
find extensive background on self-similar groups in the monograph of Nekra- 
shevych [17J, some interesting calculations particularly relevant to the basilica 
group in [ID] , and a review of the most recent developments in [18J. 

Since local regular Dirichlet forms and their Laplacians are in one-to-one 
correspondence, up to a natural equivalence, with symmetric continuous dif- 
fusion processes and their generators, our analysis allows the construction of 
diffusion processes on the basilica Julia set. Random processes of this type are 
interesting because they provide concrete examples of diffusions with nonstan- 
dard behavior, such as sub-Gaussian transition probabilities estimates. For a 
detailed study of diffusions on some finitely ramified self-similar fractals see 
[2] and references therein. The background on Dirichlet forms and Markov 
processes can be found in [6]. 

Our construction starts with providing the basilica Julia set with a finitely 
ramified cell structure, (Definition 12. II in Section [2]). According to [23], such a 
cell structure makes it possible to use abstract results of Kigami [12] (see also 
|13j ) to construct local regular Dirichlet forms that yield a topology equivalent 
to that induced from C. In Section [3] we describe all such Dirichlet forms, and 
in Section H] we describe Laplacians corresponding to these forms and arbitrary 



LAPLACIANS ON THE BASILICA JULIA SET 3 

Radon measures. Among these Laplacians, some seem more interesting than 
others. For example, there is a family of Laplacians with sufficient symmetry 
that their graph approximations admit a spectral decimation like that in [fl 
El [201 EE] ; this allows us to describe the eigenvalues and eigenfunctions fairly 
explicitly in Section [5j The energy renormalization factor is 2 and the spectral 
dimension is log 9/ log 6 in this case. 

In Section Owe describe the unique, up to a scalar multiple, Dirichlet form 
and Laplacian that are conformally invariant under the dynamical system. The 
latter Laplacian does not have spectral decimation and we cannot determine 
its eigenstructure, but its spectral dimension can be computed to be 4/3 using 
the renewal theorem in [9]. The energy renormalization factor is v^2 in the 
conformally invariant case. One can find related group-theoretic computations 
and discussions in [18] . 



2. A FINITELY RAMIFIED CELL STRUCTURE ON THE BASILICA JULIA SET 

We will construct Dirichlet forms and Laplacians on the basilica Julia set 
as limits of the corresponding objects on a sequence of approximating graphs. 
In order that we may later compare the natural topology associated with 
the Dirichlet form with the induced topology from C, we will require some 
structure on these approximations. The ideas we need are from [121 123], in 
particular the following definition is almost identical to that in [23]. The only 
change is that we will not need the existence of harmonic coordinates, and 
therefore do not need to assume that each V a has at least two elements. 

Definition 2.1. A finitely ramified set F is a compact metric space with a 
cell structure 3 = {F a } aGA and a boundary (vertex) structure V = {V a } ae ji 
such that the following conditions hold. 

(FRCS1) A is a countable index set; 

(FRCS2) each F a is a distinct compact connected subset of F\ 
(FRCS3) each V a is a finite subset of F a ; 
(FRCS4) if F a = U- =1 F aj then V a C |jj=i ^5 
(FRCS5) there exists a filtration {A n }™ =0 such that 

(i) A n are finite subsets of A, A = {0}, and F = F; 

(ii) A n R A m = if n ^ m; 

(iii) for any a G A n there are a\, au G A n +\ such that F a = 



UL F a 



'3 



(FRCS6) F a > f]F a — V a i f] V a for any two distinct a, a' G A n ; 
(FRCS7) for any strictly decreasing infinite cell sequence F ai ~D F a2 D ... there 
exists x G F such that f] n>1 F an = {x}. 

If these conditions are satisfied, then 

(F,?,V) = (F,{F a } aeA ,{V a } aeA ) 



is called a finitely ramified cell structure. 
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Notation 2.2. We denote V n = [j aeAn V a . Note that V n C V n+1 for all n ^ 
by Definition 12.11 We say that F a is an n-cell if a G A n . 

In this definition the vertex boundary Vq of F = F can be arbitrary, and in 
general may have no relation with the topological structure of F . However the 
cell structure is intimately connected to the topology, as the following result 
shows. 

Proposition 2.3 Q23J). The following are true of a finitely ramified cell struct- 
ure. 

(1) For any x G F there is a strictly decreasing infinite sequence of cells 
satisfying condition (FRCS7) of the definition. The diameter of cells 
in any such sequence tend to zero. 

(2) The topological boundary of F a is contained in V a for any a G A. 

(3) The set V* = [J aEA V a is countably infinite, and F is uncountable. 

(4) For any distinct x,y G F there is n(x, y) such that if m ^ n(x, y) then 
any m-cell can not contain both x and y. 

(5) For any x G F and n ^ 0, let U n (x) denote the union of all n-cells that 
contain x. Then the collection of open sets U = {U n (x)°} x< zF,n^o is a 
fundamental sequence of neighborhoods. Here B° denotes the topologi- 
cal interior of a set B. Moreover, for any x G F and open neighborhood 
U of x there exist y G V* and n such that x G U n (x) C U n (y) C U. In 
particular, the smaller collection of open sets W = {U n (x)°} x£ v t ,n^o is 
a countable fundamental sequence of neighborhoods. 

In general a finitely ramified fractal may have many nitrations, and the 
Dirichlet forms, resistance forms and energy measures we will discuss later are 
independent of the filtration. However it is natural in the context of a self- 
similar set to consider a filtration that is adapted to the self- similarity. We 
now define a finitely ramified cell structure and a filtration, which have certain 
self-similarity properties, on the basilica Julia set. 

By definition, the 0-cell is the basilica Julia set fractal, which we denote by 
J. Let us write a = l ~2 f° r one °f the fixed points of z 2 — 1. The interiors of 
four 1-cells are obtained by removing the points ±a; this disconnects the part 
of J surrounding the basin around into symmetric upper and lower pieces, 
and separates these from two symmetric arms, one on the left and one on the 
right, see Figure [2] (and also Figure E]). The top and bottom cells we denote 
Jm and J( 2 ) respectively, and the left and right cells we denote J( 3 ) and 
respectively. The cells J(i) and J(2) each have two boundary points, while J(3) 
and J(4) each have one boundary point. In the notation of Definition 12. 1[ 

V(i) = V(2) = {±«}, V( 3 ) = {-a}, 1/(4) = {a} 

and therefore the boundary set of the fractal is Vo = {—a, a}. Note that 
the other fixed point, b = 1+ 2 , does not play any role in defining the cell 
structure. 
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Figure 2. The the Julia set of z 2 — 1 with the repulsive fixed 
points a = 1 ~ 2 V ^ and b = 1+ 2 V ^ circled. 




Figure 3. An arc-type cell 




Figure 4. A loop- type cell 

For n > 1 we set A n = {1, 2, 3, 4} x {1, 2, 3}" -1 . To define the smaller cells, 
we introduce the following definition. If a cell has two boundary points, it is 
called an arc-type cell. If a cell has one boundary point, it is called a loop-type 
cell. 

Each arc-type n-cell J a is a union of three n + 1-cells J a i, J a 2 and J a z] J a i 
and J a2 are arc-type cells connected at a middle point, while J a3 is a loop- type 
cell attached at the same point (Figure [3]). 

Each loop-type n-cell J a is a union of three n + 1-cells, J a \, J q2 and J a3 ; 
J a i and J a 2 are arc-type cells connected at two points, one of which is the 
unique boundary point v a G V( a ), while the other is the boundary point of the 
loop-type cell J a3 (Figure 0j). 

The existence of this decomposition is a consequence of known results on 
the topology of quadratic Julia sets. In essence we have used the fact that the 
filled Julia set is the closure of the union of countably many closed topological 
discs, and that the intersections of these discs are points that are dense in J 
and pre-periodic for the dynamics. The Julia set itself consists of the closure of 
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the union of the boundaries of these topological discs. This structure occurs 
for the Julia set of every quadratic polynomial z 2 + c for which c is in the 
interior of a hyperbolic component of the Mandelbrot set or is the intersection 
point of two hyperbolic components, so in particular for the basilica Julia set 
because c = — 1 lies within the period 2 component. Details may be found 
in [HI HI HE]. These general results imply that a FRCS may be obtained for 
all quadratic Julia sets with suitable c values in the manner similar to that 
described above, however the basilica Julia set is a sufficiently simple case that 
the reader may prefer to verify directly that the existence (see, for instance, 
[161 Theorem 18.11]) of internal and external rays landing at a implies that 
deletion of ±a decomposes J into the four components Ju\, % = 1,2, 3, 4, while 
the remainder of the decomposition follows by examining the inverse images 
of these sets under the dynamics. 

Definition 2.4. The basilica self-similar sequence of graphs G n have vertices 
V n as previously described. There is one edge for each pair of vertices joined 
by an arc-type n-cell, as well as one loop at each vertex at which there is a 
loop-type n-cell. The result is shown in Figure [51 and we emphasize that it is 
highly dependent on our choice of filtration. 




Figure 5. Basilica self-similar sequence of graphs: graphs Go 
and G\. 

The sequence of graphs in Figure [5] is well adapted to the construction of the 
Kigami resistance forms, and hence the Dirichlet forms, on J. For this reason 
it plays a prominent role in Section [21 In Section a spectral decimation 
method for this sequence of graphs is used to obtain a full description of the 
corresponding Laplacian. 

It should be noted, however, that this is not the only sequence of graphs that 
we will consider. A sequence that is arguably more natural, is the conformally 
invariant graph-directed sequence of graphs for the basilica Julia set, shown 
in Figure [61 

These graphs will be considered in Section El where detailed definitions 
can be found. Their construction is related to group-theoretic results [T71 
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Figure 6. Basilica conformally invariant graph-directed se- 
quence of graphs. 



A 



B 



B 



I 

Figure 7. The substitution scheme for the basilica conformally 
invariant graph-directed sequence of graphs. 

and references therein], and in particular to the substitution scheme in 
Figure [7J The cell structure and the filtration could be defined starting with 
the single point boundary set {a}, and then taking the inverse images P~ n {a} 
of this point under the polynomial P(z), which is of course different from V n 
in Definition 12.41 More precisely, for any n and k we have V n ^ P~ h {a}, even 

though V. = U„> Vn = U>0 P~ n {*}- ' 

3. KlGAMl'S RESISTANCE FORMS FORMS ON THE BASILICA JULIA SET AND 

THE LOCAL RESISTANCE METRIC 

One way of constructing Dirichlet forms on a fractal is to take limits of resis- 
tance forms on an approximating sequence of graphs. We recall the definition 
from [11] , as well as the principal results we will require. 

Definition 3.1. A pair (£,Dom£) is called a resistance form on a countable 
set K if it satisfies the following conditions. 

(RF1) Dom£ is a linear subspace of £(V*) containing constants, £ is a non- 
negative symmetric quadratic form on Dom£, and E(u,u) = if and 
only if u is constant on V*. 

(RF2) Let ~ be the equivalence relation on Dom£ defined by u ~ v if and 
only if u — v is constant on V*. Then (£/~, Dom£) is a Hilbert space. 
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(RF3) For any finite subset V C V* and for any v G £{V) there exists u G 

Dom£ such that u\ y = v. 
(RF4) For any x, y G K the resistance between x and y is defined to be 

v( \ J (u(x)-u(y)) 2 wnl/ 

K{x, y) = sup < —, r : u G Dome, c{u, u) > (J > < oo. 

[ £{u,u) J 

(RF5) For any u G Dom£ we have the E(u,u) ^ £(u,u), where 

1 if u{x) ^ 1, 
//(.,•) = ^ u(x) if < u(x) < 1, 
if u{x) ^ 1. 

Property (RF5) is called the Markov property. 

Proposition 3.2 (Kigami, [12]). Resistance forms have the following proper- 
ties. 

(1) The effective resistance R is a metric on V*. Any function in Dom£ 
is R- continuous; in particular, ifQ is the R-completion ofV* then any 
u G Dom£ has a unique R-continuous extension to Q. 

(2) For any finite subset U C V*, a finite dimensional Dirichlet form &u 
on U may be defined by 

£u(f,f) = inf{£(fiS50 : 9 e Domfi,^ = /}. 

There is a unique g at which the infimum is achieved. The form £V 
is called the trace of £ on U, and may be written 8.jj = Trace [7 (£). If 
U\ C U 2 then = Trace^fit/J. 

Our description of the Dirichlet forms on the basilica Julia set relies on the 
following theorems. 

Theorem 3.3 (Kigami, [12J). Suppose that V n are finite subsets ofV* and that 
U^o ^« ^ s R-dense in V*. Then 

£(/,/)= lim ZvMf) 

n— >oo 

for any f G Dom£, where the limit is actually non- decreasing. Is particular, 
£ is uniquely defined by the sequence of its finite dimensional traces £y n on 

Vn- 

Theorem 3.4 (Kigami, [12]) • Suppose that V n are finite sets, for each n there 
is a resistance form £y n on V n , and this sequence of finite dimensional forms 
is compatible in the sense that each £y„ is the trace of £y n+1 on V n , where 
n = 0,1,2,.... Then there exists a resistance form £ on V* = IJ^Lo such 
that 

£(/,/) = lim e. Vn (f,f) 

for any f G Dom£, and the limit is actually non- decreasing. 
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For convenience we will write £ n (f, f) = £y„ (/,/)■ A function is called 
harmonic if it minimizes the energy for the given set of boundary values, so 
a harmonic function is uniquely defined by its restriction to Vq. It is shown 
in P2J that any function h on Vq has a unique continuation to a harmonic 
function h, and £(h,h) = E n (h,h) for all n. This latter is also a sufficient 
condition: if g G Dom£ then £ (g,g) ^ E.(g,g) with equality precisely when 
g is harmonic. 

For any function fonV n there is a unique energy minimizer h among those 
functions equal to / on V n . Such energy minimizers are called n-harmonic 
functions. As with harmonic functions, for any function g G Dom£ we have 
£n(g,g) ^ &(g,g), an d h is n-harmonic if and only if £ n (/i, h) = E(h, h). 

It is proved in [23] that if all n-harmonic functions are continuous in the 
topology of F then any F-continuous function is i?-continuous and any R- 
Cauchy sequence converges in the topology of F. In such a case there is also a 
continuous injection 9 : Q — > F which is the identity on V*, so we can identify 
Q with the the i?-closure of V* in F. In a sense, Q is the set where the Dirichlet 
form £ "lives". 

Theorem 3.5 ( [T2J 23J ) . Suppose that all n-harmonic functions are continu- 
ous. Then £ is a local regular Dirichlet form on L 2 (Q, 7), where 7 is any finite 
Borel measure on (F, R) with the property that all nonempty open sets have 
positive measure. 

Proof. The regularity of £ follows from [T2J Theorem 8.10], and its locality 
from [231 Theorem 3]. Note that, according to [121 Theorem 8.10], in general 
for resistance forms one can consider a-finite Radon measures 7. However for 
a compact set a Radon measure must be finite. □ 

The trace of £ to the finite set V n may be written in the form 

(3.1) £»(/>/)= E ^"'(/W-ZM) 2 , 

from which we define the resistance across J a to be the value r a . Note that it 
is not the same as R(v a i,v a 2). 

The values r a may be used to define a geodesic metric that is comparable 
to the resistance metric. A path from x to y in Q consists of a doubly infinite 
sequence of vertices {^q, }^ and arc-type cells Jj connecting v aj to v a . +1 , with 
lim^.oo v a . = x and lim^oo v a . = y, the limit being in the -R-topology. The 
length of the path is the sum of the resistances of the constituent cells. If x 
(or y) is in V* we permit that the sequence begins with infinite repetition of 
x (respectively ends with repetition of y) which are considered connected by 
the null cell of resistance zero, but otherwise the v a . are distinct. Let S(x,y) 
denote the infimum of the lengths of paths from x to y\ it is easy to see from 
the finitely ramified cell structure that there is a geodesic path that has length 
S(x,y). 
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Definition 3.6. We call the geodesic metric S(x, y) the local resistance metric. 
Lemma 3.7. For x and y in Q, 

^S{x,y) < R{x,y) < S(x,y). 

Proof. First consider the special case in which x and y are both in a loop- 
type cell J a , and neither is contained in any smaller loop-type cell. In this 
case none of the smaller loop-type cells affects R(x,y) or S(x,y), so we may 
replace each such loop by its boundary vertex. The result is to reduce the 
loop-type cell to a topological circle. Deleting x and y from this circle leaves 
two resistors, one with resistance S(x, y) and the other with resistance at 
least S(x,y). The resistance R(x,y) is the parallel sum of these, so satisfies 
hS(x, y) < R(x,y) < S(x,y). This case also applies if both x and y are in 
J(i) U J( 2 ) and neither is in any loop-type cell. 

To complete the proof we show that the resistance from x to y decomposes 
as a series of loop-type cells of the above form. Consider the (possibly empty) 
collection of loop-type cells that contain x but not y, and order them by 
inclusion, beginning with the largest. Let v ao , v ai , v a2 . . . be the boundary 
vertices of these loops, and observe that v a . — > x. Do the same for the loop- 
type cells containing y but not x, labeling the vertices f a _i, v a _ 2 , . . . . If x 
is in V* then the sequence will terminate with an infinite repetition of x, and 
similarly for y. Notice that deleting any of the v a . disconnects v aj _ 1 from v aj+1 . 
This implies both that that the effective resistances R(v aj ,v aj+l ) sum in series 
to give the effective resistance R(x,y), and that the resistances S(v aj , v aj+1 ) 
sum to S(x, y). However each of the configurations v a ,, v a , +1 is of the form of 
the special case given above, so satisfies 

-S(v aj ,v aj+1 ) < R(v a] ,v aj+1 ) < S{v aj ,v aj+1 ). 
Summing over j then gives the desired inequality. □ 

By virtue of Theorem 13.31 and ( 13. ip it is apparent that we may describe a 
resistance form on V* in terms of the values r a . The simple structure of the 
graphs makes it easy to describe the choices of {r a } a that give a resistance 
form. 

Lemma 3.8. Defining resistance forms on each V n by (13. ip produces a se- 
quence £ n that is compatible in the sense of Theorem \ 3.4\ if and only if for 
each arc-type cell J a , 

r a = r a i + r a2 . 

Proof. Resistance forms satisfy the well-known Kirchoff laws from electrical 
network theory (see [IT] . Section 2.1). If J a is an arc-type cell then J a \ and 
J a2 connect v a i and v a2 in series. The resistance in V\ a+ i\ between t> al and v a2 
neglecting J \ J a is then r a \ + r a2 , so is compatible with the resistance across 
J a in V\ a \ if and only if r a = r al + r a2 . In the alternative circumstance where 
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J a is a loop-type cell, there is only one boundary vertex, so r a is not denned 
and no constraint on r al and r a2 is necessary. □ 

According to Lemma 13.81 and Theorem 13.41 one may construct a resistance 
form on V* simply by choosing appropriate values r a . It is helpful to think of 
choosing these values inductively, so that at the n-th stage one has the values 
r a with | a | = n. In this case there are two types of operation involved in 
passing to the (n+ l)-th stage. For arc-type cells J a with |a| = n one chooses 
r a ± and r a2 so they sum to r a . For loop-type cells J a one chooses r a i and r a2 
freely. 

This method provides a resistance form on V* and its .R-completion Q, but 
our goal is to describe Dirichlet forms on the fractal J. In order that Q = J, 
or equivalently that the topology from C coincides with the i?-topology on V*, 
we must further restrict the values of r a . In the theorem below, S — Diam (O) 
denotes the diameter of a set O with respect to the local resistance metric 
S(x,y). 

Theorem 3.9. The local regular resistance forms on V* for which Q = J 
and the R-topology is the same as the induced C-topology are in one-to-one 
correspondence with the families of positive numbers r a , one for each arc-type 
cell J a , that satisfy the conditions 

(3.2) r a = r al + r a2 

(3.3) lim max S - Diam (J a ) = 0. 

A sufficient but not necessary condition that implies (13.31) . and is often more 
convenient, is 

(3.4) > maxr„ < cx3. 

n 

Proof. Lemma 13.81 shows that the condition (13. 2p on the r a is equivalent to 
compatibility of the sequence of resistance forms, which is necessary and suf- 
ficient to obtain a resistance form on Q by Theorems 13.31 and 13.41 

Recall that K, is C-dense in the complete metric space J, so J is the C- 
completion of V*. Similarly, Q is by definition the i?-completion of V*. Then 
Q = J and the i?-topology is the same as the induced C-topology if and only 
if every C-Cauchy sequence in V* is i?-Cauchy and vice-versa. 

Suppose there is an arc-type cell J a with r a = 0. Then the sequence defined 
by x 2 j = v a i and x 2 j+\ = v a2 is -R-Cauchy but not C-Cauchy. Conversely 
suppose there is a sequence that is -R-Cauchy but not C-Cauchy. Compactness 
of J allows us to select two distinct C-limit points x and y and these will have 
R(x,y) = 0. Then S(x,y) = by Lemma [3 .7\ thus there is a non-trivial path 
joining x to y such that r a . = for all arc-type cells J aj on the path. It follows 
that .R-Cauchy sequences are C-Cauchy if and only if r a > for all arc-type 
cells J n . 
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An equivalence class of C-Cauchy sequences is a point x G J, and as noted in 
Proposition ^. 31 x is canonically associated to the nested sequence {U n (x)j , 
where U n (x) is the union of the n-cells containing x. Hence C-Cauchy se- 
quences are -R-Cauchy if and only if for each x the resistance diameter of 
U n (x) goes to zero when n —>■ oo. Clearly this is true if 

(3.5) lim max LR-Diam (J a ) = 0. 

n-+oc a<=A n V / 

and conversely if (13.51) fails then there is e > and for each n a cell J a with 
| at | = n and .R-diameter at least e, so compactness of J gives a C-limit point 
x at which the i?-diameter of U n (x) is bounded below by e independent of n. 
Applying Lemma 13.71 we then see that C-Cauchy sequences are i?-Cauchy if 
and only if (I3.3P holds. 

For any cell J a and any x G J a there is a path from v a \ to x which contains 
at most one arc- type cell of each scale less than |a|, so the condition (13.41) 
implies f !3.3j) . To see this condition is not necessary we consider {r a } as follows. 
Fix a collection of loop-type cells J a , one for each scale |a| > 1, with the 
property that if J a is in this collection then no loop-type ancestor of J a is in 
the collection. For example the cells with addresses (3), (13), (113), (1113), .... 
If J a is in this collection set r al = r a2 — If J Q is a loop-type cell not in 

this collection set r a \ = r a2 = 2~l a L Also let r al = r a2 = r a /2 if J a is arc-type. 
For these values r a we see that any local resistance path contains at most one 
arc-type cell from this collection, so the S'-diameter of a cell of scale n > 2 is 
at most (n — 1) _1 + 2~ n ~ 2 and (13. 3p holds. However max agyln r Q = (n — 1) _1 
for each n > 2, so (13.41) fails. □ 

Corollary 3.10. Under the conditions of Theorem all the functions in 
dom(£) are continuous in the topology from C. 

Proof. It follows from (RF4) in Definition 13.11 that functions in dom(£) are 
^-Holder continuous in the i?-topology. □ 

The n-harmonic functions have a particularly simple form when written with 
respect to the local resistance metric. 

Theorem 3.11. An n-harmonic function is piecewise linear in the local resis- 
tance metric. 

Proof. The complement of V n is the finite union of cells J a with a G A n . If / is 
prescribed at v a i and v a2 then its linear extension to J a in the local resistance 
metric has f{v a3 ) satisfying r a f(v a3 ) = r al f(v a2 ) +r a2 f(v al ). However the 
terms in the trace of £ to V n +i that correspond to J a are 

^ai(/Ki) - f(v*3)) 2 + r~ 2 (f(v a2 ) - f{v a3 )) 2 

and it is clear this is minimized at precisely the given choice of f(v a3 ). We have 
therefore verified that an n-harmonic function extends from V n to V n+ \ linearly 
in the local resistance metric, and the full result follows by induction. □ 
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It is sometimes helpful to think of the local resistance metric as correspond- 
ing to a local resistance measure u, which is defined as follows. 

Definition 3.12. The local resistance measure of a compact set E is given by 



One can easily see that u has a unique extension to a positive, possibly 
infinite, Borel measure, and that S(x, y) is the smallest measure of a path 
from x to y. 

This measure has a natural connection to the energy measures corresponding 
to functions in dom(£). If £ is local and / G dom(£), then the standard way 
to define the energy measure vj is by the formula 



for any bounded quasi-continuous g G dom(£), see for instance [6]. If E is 
open, then another way to define Vf{E) is to take the limit defining £ from the 
resistance form as in Theorem 13.41 and restrict to edges in E. One may infor- 
mally think of the energy measure Vf{E) of a set E as being Vf(E) = £(/e), 
where is equal to u on E and zero elsewhere, though this intuition is non- 
rigorous because f& may fail to be in the domain of £. If h m is the piecewise 
harmonic function equal to u on V m then Uh m — > ff, and Theorem 13.111 ensures 
h m has constant density -4j on every sufficiently small cell. This sequence of 
piecewise constant densities is bounded by £(/) in L 1 (<iz/), and is a uniformly 
integrable submartingale. The limit is the density of Vf with respect to u, 
hence all energy measures are absolutely continuous with respect to v. If we 
let {fl,} be the bounded Fatou components of the polynomial P(z) = z 2 — 1 
and note that S(x, y) provides a local parametrization of the topological circle 
dQj, then the above argument gives the following description of the resistance 
form and the energy measures. 

Theorem 3.13. Under the conditions of Theorem \ 3.9[ the domain dom(£) 
of £ consists of all continuous functions such their restriction to each dQj is 
absolutely continuous with respect to the parametrization by the local resistance 
metric, and the naturally defined derivative % is square integrable with respect 
to v. Moreover, each measure v$ is absolutely continuous with respect to v, 



(3.6) 




j 



j 




dvf = ( df\ 2 
dv \dSJ 



v-almost everywhere, and 
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Note that the derivative % can be defined only up to orientation of the 
boundary components dQj, but the densities and integrals in this theorem are 
independent of this orientation. In general, v is non-atomic, cx-finite, and for 
each j we have 

< is(dilj) < oo. 

It should also be noted that v plays only an auxiliary role in this theory, and 
is not essential for the definitions of the energy or the Laplacian. 

Two specific choices of v corresponding to resistance forms of the type de- 
scribed in Theorem 13.91 will be examined in more detail in Sections and El 
In both these cases v is not finite, but cx-finite. 



4. Laplacians on the basilica Julia set 

As is usual in analysis on fractals, we use the Dirichlet form to define a 
weak Laplacian. If /j, is a finite Borel measure on J, then the Laplacian with 
boundary behavior B is defined by 

(4.1) E(f,g) = - J{A B f)gd(i for all 3 G dom B (£) 

where doing (£) is the subspace of functions in dom(£) satisfying the boundary 
condition B. In particular, if there is no boundary condition we have the 
Neumann Laplacian A^v and if the boundary condition is that g = on Vo we 
obtain the Dirichlet Laplacian Ad. We may then define a boundary operator 
such that gU) 

can be extended to a general Gauss-Green formula. 



(4.2) E(f,g) = - j (A B f)gdpL+J29(x) d nf(x) for all 9 G dom(£). 

Proofs of the preceding statements may be found in fF2\ . 

The Laplacian may also be realized as a renormalized limit of Laplacians 
on the graphs G n by using the method from [TT] . For x G V n let ip™ denote 
the unique n-harmonic function with ip^iv) = ^,2/ fo r V e where 5 x ^ y is 
Kronecker's delta. Since this function is n-harmonic, 8,(u,ip™) = £ n (u,ijj™) for 
all u G dom(£). From this and (13.11) we see that if x is in V n -\ then 

£n(«,C) = r xy{ U ( X ) - U (V)), 

y~ n x 

where y ~ n x indicates that y and x are endpoints of a common arc-type 
n-cell, and r xy is the resistance of that cell. We may view the expression on 
the right as giving the value of a Laplacian on G n at the point x 

(4.3) Alu{x)=Y,r- y l {u{x)-u{y)) 
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where the superscript r in A r n indicates its dependence on the resistance form. 
By the Gauss-Green formula 14.21 



so that 

(4.4) ( / c = ~ I( f A ;Jt; d " - - A «<*> 

as ri — > oo, which expresses A as a limit of the graph Laplacians A^, renor- 
malized by the measure \l. 



5. Spectral decimation for a self-similar but not conformally 

invariant laplacian 

The procedure in ( 14. 3f) and (14.41) is especially of interest when both the 
resistance form and the measure have a self-similar scaling that permits us to 
express A^ in terms of the usual graph Laplacian 

A n u(x) = ^2 ( u ( x ) - u (y)) 

and to simplify the expression for the measure. Consider for example the 
simplest situation, in which a resistance form is constructed on J by setting 
r a = 2~l a l, where \a\ is the length of the word a and using (13.11) . and a Dirichlet 
form is obtained as in Theorem 13.41 We take the measure /zb to be the natural 
Bernoulli one in which each n-cell has measure (4 ■ S^ 1 )^ 1 for n > 1. In this 
case (14.31) simplifies to A r n u(x) = 2 n A n u(x), and since J ip^dfiB = 2~ 1 3~ n we 
may reduce (I4.4p to 

(5.1) 2 • 6 n A n u(x) -> -Au{x). 

The negative sign occurring on the right of (15. ip is a consequence of the fact 
that A n is positive definite, whereas the definition (14.11) gives a negative def- 
inite Laplacian. The former is more standard on graphs and the latter on 
fractals. 

For the remainder of this section we study the particular Laplacian defined in 
(15. ip using its graph approximations. We begin by computing the eigenstruc- 
ture of the graph Laplacian on G n using the method of spectral decimation 
(originally from [5J [201 EI], though we follow [U US])- The situation may be 
described as follows. The transition matrix M n for a simple random walk on 
G n is an operator on the space of functions on V n . If we decompose this space 
into the direct sum of the functions on V n -\ and its orthogonal complement, 
then M n has a corresponding block form 



(5.2) 



16 



LUKE G. ROGERS AND ALEXANDER TEPLYAEV 



in which the matrix A n is a self-map of the space of functions on V n -i. De- 
fine the Schur complement S to be A n — B n D~ l C n , and consider the Schur 
complement of the matrix M n — z = M n — zl: 



(5.3) S n (z) = A n -z- B n (D n - z)- x C n . 
If it is possible to solve 

(5.4) S n (z) = <j )n (z)(M n _ 1 -R n (z 

where 4> n {z) and R n (z) are scalar- valued rather than matrix- valued rational 
functions, then we say that M n and M n _i are spectrally similar. If we have 
a sequence {M n } in which each M n is a probabilistic graph Laplacian on G n 
and M n is spectrally similar to M n _i, then it is possible to compute both 
the eigenvalues and eigenfunctions of the matrices M n from (j) n (z) and R n (z). 
Excluding the exceptional set, which consists of the eigenvalues of D n and the 
poles of (pn(z), it may be shown that z is an eigenvalue of M n if and only if 
Rn{z) is an eigenvalue of M n _i, and the map / i— > / — (D n — z)~ 1 C n f takes 
the eigenspace of M n _i corresponding to R n (z) bijectively to the eigenspace 
of M n corresponding to z ( |15|, Theorem 3.6]). 

Now consider a self-similar random walk on the graph G n in which the 
transition probability from u a3 to v a \ in an arc-type cell is a fixed number 
p G (0, 1/2). The transition matrix for the cell J a has the form 



(5.5) M 



The results of pQ imply that the spectral decimation method is applicable to 
the graph G n . Moreover, self-similarity implies that both <j) n (z) and R n {z) are 
independent of n and may be calculated by examining a single cell J a . From 
(15. 5p we see that the eigenfunction extension map is 

(D-z)- 1 C={ 3 fr s 

meaning that the value at v a3 of a A| Q | +1 eigenfunction is -^—^ times the sum 
of the values at v a \ and v a 2- Since M n _i on a single cell is simply 




M 

we find that 



1 -1 
-1 1 



(j)(z) = — - — , and 
2p — z 

(5.6) R(z) = ^^z - -z\ 

p p 
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The exceptional set is exactly the point {2p}. If we choose the initial Laplacian 
on G to be 

A =f q ~ q ) 

for some < q < 1 then we can apply Proposition 4.1 and Theorem 4 of [1] 
to compute both the multiplicities and the eigenprojectors. 

Theorem 5.1. The eigenvalues of the Laplacian A n on G n are given by 

a(A ) = {0,2g}, 

a{K) = ( U R ' m {M) U (^{O' 2 *?}) • 

\m=0 / 

Moreover, if z G R~ n {0,2q} then mult n (z) = 1 and the corresponding eigen- 
functions have support equal to J; if z G R~ m {2p} then mu\t n (z) = 2 ■ ^ n ~ m ~ 1 
and the corresponding eigenf unctions vanish on V n - m -\- 

Proof. For z G R~ n {0, 2q} the result follows from Proposition 4.1 (i) and The- 
orem 4(i) of pp. For z G R~ m {2p} the result follows from Proposition 4.1(iii) 
and Theorem 4(iii) of [lj. In particular, 

mult„(2p) = 4 • - |K-i| + mult n _i(i2(2p)), 
where R{2p) = 2, which is not in the spectrum of any A^, and \Vk\ = 2-3 k . □ 

Corollary 5.2. The normalized limiting distribution of eigenvalues (also called 
the integrated density of states ) is a pure point measure k with atoms at each 
point of the set 

oo 

U R' m {M, 

m=0 

Moreover, if z G R~ m {2p] then k({z}) = 2-3~ m ~ 1 . There is one atom in each 
gap of the Julia set of R. 

A special case occurs if we make the convention that every edge can be 
traveled in both directions with equal probability, in which case each of the 
G n is a regular graph of degree 4. This simple random walk has p = q = \ 
from which R(z) = 6z — 4z 2 . Since our graphs have 2 • 3 n vertices, we conclude 
that in this case the spectral dimension of the corresponding infinite graphs is 

2 log 3 

d s = . 

logo 

One may also consider weighted Laplacians on the infinite graphs by varying 
the parameter p. 

We saw at the end of Section H] that the Laplacian A on the fractal J may be 
obtained as a limit of graph Laplacians A n , provided that both the Dirichlet 
form and the measure have self-similar scaling. Under these circumstances, 
the spectral decimation method gives a natural algorithm for constructing 
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eigenfunctions of the Laplacian on the fractal. This method was first developed 
for the Sierpinski Gasket fractal [201 ED E] . 

We illustrate this method for the special self-similar case where the resis- 
tance form on J satisfies (13.1 j) with 

_ o-H 

'a — " i 

where \a\ is the length of the word a, and the Dirichlet form is obtained using 
Theorem 13.41 In this case 

where A n is the graph Laplacian on G n with equal weight \ on each edge. 
This is equivalent to setting p = q = |. Correcting for the extra factor of | in 
the graph Laplacian we find from (15.11) 

(5.7) 8 • 6 n A n u(x) -> -Au(x). 

Here we take that the measure \x in (14. ip is the the natural Bernoulli measure 
Hb for which each n-cell has /^-measure equal to (4 ■ 3" -1 )- 1 for n > 1. 

Now suppose that {u n } is a sequence of eigenfunctions of A n with eigen- 
values A n , and the property that u n = u m on V m for m < n. Further assume 
that 6 n A n converges and that the function u defined on V* by u(x) = u n (x) for 
x G V n is uniformly continuous, and thus can be extended continuously to J. 
Then (15 .7p implies that u is a Laplacian eigenfunction on J with eigenvalue 
A = — 81im6 n A n . From the formula ( 15. 6p for R we have 

3 + e n> /9-4A n _i 
An 

where e n is one of ±1 for each n. If only finitely many e n equal +1 then 6 n \ n 
converges and it is easily verified that u is uniformly continuous on V*, so this 
method constructs a large number of eigenfunctions. It is actually the case 
that it constructs all eigenfunctions, though we will only show this for the 
Dirichlet Laplacian. 

The Dirichlet eigenfunctions corresponding to R~ m (2p) = R~ m (^) produce 
Dirichlet eigenfunctions on J. Via an argument from [5] , this provides a precise 
description of the spectrum of the Dirichlet Laplacian A^. Let ijj(x) = 3 ~^^~ 4x 
and 

*(x) = lim 6> n (x) 

n— >oo 

in which the limit is well-defined on a neighborhood of zero by the Koenig's 
linearization theorem (see [16]). Note that ^(0) = and ^'(0) = 1, so that ^ 
is also invertible on a neighborhood of 0. In the above construction of Dirichlet 
eigenvalues we asked that all but finitely many of the inverse branches of R 
be exactly ip, so that for any such A = — 81im6 n A n there is n such that 
A n+ i = ip(X n ) for all n > uq. It follows that 

A = -8 ■ lim 6 no 6 n - no V n - no (A no ) = -8 ■ 6 n °*(A no ) 
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where A no = R m (|) for some < m < n . 

Theorem 5.3. The spectrum of Ad on J consists of isolated eigenvalues 

A = -8-6 n °#(iT m (~)) 

with multiplicity 2-3 n °~ m ~ 1 , for each no > 1 and < m < uq. The correspond- 
ing eigenf unctions are those obtained from the eigenfunctions in Theorem \5.1\ 
by spectral decimation. 

Proof. Kigami [12] proves that there is a Green's operator with a kernel that is 
uniformly Lipschitz in the resistance metric, hence the resolvent of the Lapla- 
cian is compact and the spectrum of the Laplacian is discrete (pure point with 
isolated eigenvalues of finite multiplicity accumulating to infinity). Since A D is 
negative definite, the spectrum consists of a decreasing sequence Xj of negative 
real eigenvalues that accumulate only at — oo. 

We have seen that the spectral decimation construction produces some 
Dirichlet eigenvalues and their eigenfunctions. The standard way to determine 
that all points in the spectrum arise in this manner is a counting argument 
due to Fukushima and Shima [3]. As the argument holds essentially without 
alteration, we only sketch the details. 

Expanding the Green's kernel g(x, y) of the Laplacian as an L 2 -series in the 
eigenfunctions, we find that 

- J g(x, x)d/i B (x) = J; 

i 

where the sum is over the eigenvalues of Ad, each repeated according to its 
multiplicity. Similarly, if we let g m be the Green's kernel for —8 ■ Q m A m and 
let \x m be the measure with equal mass at each point of V m , then 

— / 9m(%i x)d^L m {x) = ^ ^ 
J 3 K 3 

where the sum is over its eigenvalues. However g is continuous and equal to 
g m onV m , and the measures /i m converge weak* to /ie, so as m — > oo the sum 
of all -jijy converges to the sum of j-. 

Now each is —8 • 6 m \ ( j n \ where G i?~" m (|), and any sequence 

—8 • Q m X ( j n ^ satisfying the conditions of the spectral decimation algorithm 
converges to some eigenvalue Aj of Ajj. With a little care it is possible to show 
that Y\. -4ry converges to the sum V, y~, over those eigenvalues that arise 

from the spectral decimation. Since Y\- -7=7 also converges to -r-, we conclude 
that the spectral decimation produces all eigenvalues. □ 

It is worth noting that eigenfunctions also have a self-similar scaling prop- 
erty. Specifically, let f a denote the natural map from J(i) to J a if J a is an 
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arc-type cell, and from J( 3 ) to J a if J a is a loop-type cell. This natural map is 
defined in the obvious way on the boundary points and then inductively ex- 
tended to map V n H J(i) to H J a (respectively V n fl J( 3 ) to V^ + | a |_i fl J a ) 
for each n, whereupon it is defined on the entire cell by continuity. By the 
definition of the Dirichlet form, this composition scales energy by 2 1 -!"!, and 
by ( 15.11) it scales the Laplacian by 6 1_ ' a '. More precisely, if u is such that 
(A — X)u = then A(u o f a ) = 6 1_ ' a '(An) o / Q , so A(u o f a ) is a Laplacian 
eigenfunction with eigenvalue 6 1_ ' a 'A. 

The scaling property provides a very simple description of the Dirichlet 
eigenfunctions. Suppose u is a Dirichlet eigenfunction obtained as the limit of 
u n according to the spectral decimation, and let m be the scale with X m = \- 
Then u m vanishes on V m -x, so if |a| = m then u m o f a is a Dirichlet eigenfunc- 
tion on J(!) (or J( 3 )) with eigenvalue 6 1_m . There is a one dimensional space of 
such functions (note that whether the function is on Jm or Jt 3 \ is immaterial 
because it vanishes on the boundary), spanned by the Dirichlet eigenfunction 
on J(i) with value 1 at V13. It follows that the Dirichlet eigenfunctions are 
all built by gluing together multiples of this function on individual cells of a 
fixed scale m, subject only to the condition that the values on V m give a graph 
eigenfunction with eigenvalue |. 

6. CONFORMALLY INVARIANT RESISTANCE FORM AND LAPLACIAN 

In this section we decompose J as a union of a left and right piece J = 
JiU Jr, where 

J L = Jn{z:Re(z)^ 1 -fi} = J (3) 
J R = J n {z : Re(z) ^ ^} = J {1) U J (2) U J {4) . 

The sets meet at a = ' , which is the fixed point of P(z) = z 2 — 1. The 
polynomial P(z) maps Jl onto Jr by an one-to-one mapping, and the piece 
J(4) C Jr onto Jr by a one-to-one mapping. It also maps the central part 
J(i) U J(2) of Jr onto J L by a two-to-one mapping. Therefore the following 
directed graph 

J l ^ Jr 1 

corresponds to the action of P(z), and defines a graph directed cell structure on 
J. Note that K, = U m P _m {a} and that the preimages of arc-type cells under 
P are also arc-type cells, while the preimages of loop-type cells are loop-type 
cells except in the case of J( 3 ) = Jl for which the preimages are J(i) and J( 2 ). 
This construction is related to group-theoretic results about these graphs and 
[TTt ITS! and references therein] , and in particular to the substitution scheme in 
Figures [6] and d in which the labeling of components is Jl = A and Jr = B. 

As usual, we are interested in Dirichlet forms and measures that have a self- 
similar scaling under natural maps of the fractal. In this case, the mapping 
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properties described above show that if £ is a Dirichlet form on J then we 
may define Dirichlet forms £* on the cells Jr^, i = 1,2,3,4 by setting 

£*(«) = £(u o P) for u on Ju\ with moPg dom(£) 

where u o P is taken to be zero off P(J^) in each case. The form £ is then 
self-similar under the action of P if for u G dom(£) 

(6.1) £( M ) = p^£ i ( M | J(i) ) 

i 

for some p. 

Theorem 6.1. Among the resistance forms identified in Theorem \3.9i there is 
one that has a self-similar scaling under the action of P(z) and is symmetric 
under complex conjugation. It is unique up to a scalar multiple, and has scaling 
factor p = y/2. 

Proof. By Theorems 13.31 and 13.41 a necessary and sufficient condition for (16.11) 
to be true is that the trace of both £ and J2i to V m are equal for each m. 
The trace of £ to V m is a resistance form 

£ m («) = ^ r a l { u { v al) -u(v a2 )) 2 

as in (13.11) . The trace of Y2i ^ to * s f° un d by minimizing the energy 
when values on V m are fixed, and each £* may be minimized separately. Thus 
for each i the restriction of u to Ju\ has the property that u o P is energy 
minimizing on P(J^). The result is therefore a resistance form in which the 
resistance across an arc-type cell J a is equal to the resistance of the form £ 
across P(J/ a )). 

We conclude that ( 16. ip is true if and only if £ is the limit of resistance forms 
£ m with rp(j (ct) ) = prj,y There is only one value of p for which this can be 
satisfied. To see this, note that P 2 maps both J(u) and J( 2 2) to J(i), and both 
J ( i2) and J ( 2i) to J( 2 ), so r (11) = r (22 ) = p~ 2 r(i) and r (12) = r (2 i) = P _2 ^(2)- 
However r (U ) + r ( i 2 ) = r (i ) and r (2 i) + r (22 ) = r {2) , so r (1) = r (2 ) and p 2 = 2. 
Also r(ji) and rus) are equal to \r(i) for z = 1, 2. 

Observe that for any arc-type cell Ji a ) there is a unique m = m(a) so 
P m (J(a)) = for one of i = 1,2. According to our reasoning thus far, we 
must have rn\ = r( 2 ) and r a = 2~ m( - a ^ 2 r^iy It remains to be seen that these 
resistances satisfy the conditions of Theorem [331 The condition limiai^oo r a = 
is immediate, and one can easily verify (13.31) . in particular by the computation 
in Theorem l6.2[ For any a we have P m ( Q ) ( J( a i)) and P m ( a ) ( J( a 2)) are J(n) and 
J(i2) for one of i — 1, 2, so the second condition r a \ + r a2 = r a is equivalent to 
+ ^12 = i~x, and the latter has already been established. □ 

Recall that we have a local resistance metric and a measure dv (see Defi- 
nition 13.61 (13. 6p and Theorem 13.131) corresponding to a sequence of resistance 
forms. For the resistance form from Theorem 16.11 the measure dv is related 
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to the harmonic measure on Fatou components. In this case the measure is 
infinite. It is well known that in our situation each Fatou component DcC 
is a topological disc with locally-connected boundary, so the Riemann map a 
from the unit disc to Q has a continuous extension to the unit circle. The 
harmonic measure from the point x G Q is the image on dQ of the Lebesgue 
measure on the circle under the Riemann map with i— > x. It is the same as 
the exit probability measure for a Brownian motion in Q started at x, and is 
the representing measure for the linear functional on C(dQ) that takes / to 
the value of the harmonic extension of / at x. 

Theorem 6.2. Let dv be the measure corresponding to the unique resistance 
form from Theorem \ 6.1\ with normalization rp) = ~ . Then v is an infinite 
measure that has the self-similar scaling v(P(E)) = y/2v(E) for any set E on 
which P is injective. 

Let Q = Q be the Fatou component of P that contains the critical point 
0, and for each other bounded component Qj of the Fatou set of P let rrij be 
the unique number such that P mj maps Qj bijectively to Q. Let di/j be the 
harmonic measure on dQj from the point P~ mj (0) G ttj. Then 



Proof. Let a be a Riemann map from the unit disc to Q with cr(0) = 0. Since 
P 2 is a two-to-one map of Q onto itself, a -1 o P 2 o a is a two-to-one map of 
the unit disc onto itself. A version of the Schwartz lemma (for example, that 
in [Tj5]) implies that a~ x o P 2 o a = cz 2 , where c is a constant with |c| = 1. 
Moreover, there is a unique a such that o~(l) = a, and so o~~ 1 oP 2 oo~ = z 2 . Since 
Q is locally connected the Riemann map extends to the boundary. Pulling back 
the measure dv to the circle via o gives a Borel measure that scales by 2 under 
zhz 2 . Consider the set of 2 m preimages of a under the composition power 
P 2m that lie in dil. These preimages divide v\ dU into 2 m equal parts. The 
preimages of these 2 m points under a are binary rational points on the unit 
circle that divide the Lebesgue measure into 2 m equal parts. It follows that dv 
is a multiple of the harmonic measure for the point cx(0) = 0, and since they 
both have measure 1 they are equal. A similar alternative construction is to 
consider an "internal ray" which is the intersection of the negative real half- 
line with Q, and its preimages. Then the harmonic measure can be determined 
in the usual way by computing angles between these rays. 

The bounded Fatou components of P are Q and the topological discs en- 
closed by loop-type cells. The argument we have just given applies to any 
such component fL, except that the Riemann map is P~ m i o a, so dv\ n is 

a multiple of the harmonic measure dvj from the point P~ mj (0). The re- 
sult then follows from the proof of Theorem 16.11 where it is determined that 




j=0 



v(dn 3 ) 



2~ rrij/2 



□ 
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It is natural to compare this to other measures on J. We saw in Theo- 
rem 13.131 that the energy measures are absolutely continuous to dv. Another 
standard measure to consider is the unique balanced invariant probability mea- 
sure of P, denoted \ip. It can be obtained, for instance, as the weak limit of the 
sequence of probability measures /i m , where each \i m is 2~ m times the counting 
measure on the 2 m preimages of a. An alternative construction of \ip defines 
it as the harmonic measure from infinity, which also can be determined in the 
usual way by computing angles between the external rays. This measure is a 
Bernoulli-type measure that has the self-similar scaling fip(P(E)) = 2fip(E) 
for any set E on which P is injective (or any set if we incorporate multiplicity). 
The measures /ip and v are singular, as may be verified directly by comparing 
Hp to v. Indeed, /ip has measure 2~ m on the preimages P~ m {J^)), i = 1,2 
whereas v(P~ m ( J (i) )) = 2~ m / 2 . 

Let Ap be the Laplacian corresponding to the unique conformally invariant 
£ and the balanced invariant measure /ip. Because of Theorem 16. 1\ Ap is (up 
to a constant multiple) the only Laplacian that has self-similar scaling under 
the action of P(z) = z 2 — 1, and its scaling factor is 2v^2- 

4 

Theorem 6.3. The spectral dimension of Ap is equal to — . 

Proof. Since J has graph-directed fractal structure, the method of [HI [TJ] is 
applicable. This reduces the spectral dimension computation to finding s such 
that the spectral radius of the matrix 

(2V2) 

is equal to one. Thus s = | and d s = 2s = |. □ 
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